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The present-day water cycle on Mars has implications for habitability 

and future human exploration. Water ice clouds and water vapour have 
been detected above the Tharsis volcanic province, suggesting the active 
exchange of water between regolith and atmosphere. Here we report 
observational evidence for extensive transient morning frost deposits on 
the calderas of the Tharsis volcanoes (Olympus, Arsia and Ascraeus Montes, 
and Ceraunius Tholus) using high-resolution colour images from the Colour 
and Stereo Surface Imaging System on board the European Space Agency’s 
Trace Gas Orbiter. The transient bluish deposits appear on the caldera 

floor and rim in the morning during the colder Martian seasons but are not 
present by afternoon. The presence of water frost is supported by spectral 
observations, as well as independent imagery from the European Space 
Agency’s Mars Express orbiter. Climate model simulations further suggest 
that early-morning surface temperatures at the high altitudes of the volcano 
calderas are sufficiently low to support the daily condensation of water—but 
not CO,—frost. Given the unlikely seasonal nature of volcanic outgassing, 

we suggest the observed frost is atmospheric in origin, implying the role 

of microclimate in local frost formation and acontribution to the broader 


Mars water cycle. 


The Tharsis Rise is a large volcanic province in the tropics of Mars’ 
(latitude range: +40° N, longitude range: 220-300° E). It is a broad 
topographic dome that rises about 5 km above the surrounding 
terrain and covers a region 5,000 km wide’. It contains some of the 
Solar System’s largest and tallest volcanoes’, such as Olympus Mons 
(21 km altitude), Arsia Mons (18 km), Ascraeus Mons (18 km) and 
Pavonis Mons (14 km), but also smaller shield volcanoes such as Cerau- 
nius Tholus (9 km). Volcanic activity on Mars has been concentrated 
predominantly in this region throughout the planet’s geological his- 
tory, persisting into current times, as evidenced by lava flows that are 


as recent as 2.4 million years old*. No current volcanic activity has been 
detected in Tharsis, although recent geophysical data show that Mars 
is still geodynamically active’ ’. 

Notable orographic water ice clouds and other atmospheric 
phenomena have been observed in Tharsis®* ’. Water ice clouds play 
a fundamental role in cycling water on Mars, moving moisture for 
thousands of kilometres from polar regions to relatively dry equa- 
torial areas®*. In addition, Tharsis is situated along the route of an 
important cross-equatorial exchange of water vapour, where approxi- 
mately 10” kg of water is annually transferred between the northern 
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Fig. 1| Frost detection on the Olympus Mons caldera by CaSSIS, HRSC and 
NOMAD. a, Global view of Mars with white box marking the location of Olympus 
Mons. b, HRSC wide-angle image of Olympus Mons acquired in the early morning 
(LST = 7:20, Ls = 346.7°, latitude = 18.2° N, longitude = -133.2° E). The black 
dashed line indicates the orbit of the TGO corresponding to the images in d and 
e. The white box highlights the close up inc. c, Zoomed-in view of the Olympus 
Mons caldera. The white and blue dashed rectangles show the footprints of 

the CaSSIS and NOMAD-LNO observations, respectively. d, High-resolution 

(4.5 m pixel) CaSSIS colour image of frost on the caldera floor and northern 
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rim of Olympus Mons (LST = 7:11, Ls = 344.1°). Frost is absent on the well-lit 

steep slopes. The blue rectangle marks the footprint of the one NOMAD-LNO 
observation that falls within the frost-covered area. e, NOMAD-LNO channel 
observation of the Olympus Mons caldera. The ice index values (Methods) 
indicate the presence of frost over the caldera floor (> + 30). The coloured areas 
onthe plot indicate the confidence intervals. HRSC image ID: hn889_0000 (b,c). 
CaSSIS colour image ID: MY36_022332_162_0_NPB (d). NOMAD-LNO observation 
ID: 20221125 082524 (e). Credit: b, ESA/DLR/FU Berlin; d, ESA/TGO/CaSSIS under 
a Creative Commons license CC-BY-SA 3.0 IGO. 


and southern hemispheres through the solstitial Hadley cells». Atmos- 
pheric observations’ have revealed a localized enrichment in water 
vapour above the Tharsis volcanoes, suggesting that an active exchange 
of water vapour between the regolith and the atmosphere may be 
ongoing, probably facilitated by desorption from the regolith and/or 
sublimation of frost. A subsequent study” confirmed the water-vapour 
enrichment over these areas but hypothesized that the local circu- 
lation pattern typical of the volcanic region is possibly responsible 


for the enrichment as it may carry considerable amounts of water 
vapour upslope. 

Apart from the polar regions, water ice manifests on the sur- 
face as seasonal frost in mid- and low-latitude locations. NASA’s (the 
National Aeronautics and Space Administration’s) Viking 2 lander 
detected water frost at -48° N’S°. In addition, orbital observations 
from a variety of instruments revealed that water frost can occur up 
to 13° S in the southern hemisphere and as low as 32° N on shaded 
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Fig. 2 | Frost occurrence ona crater in the Arsia Mons caldera. a, Frost on the 
shadowed slope of the crater in an early-morning observation during southern 
winter in MY 35 (latitude = —8.74° N, longitude = -121.14° E). b-d, No frost in 

an early-morning observation (b) and no frost in afternoon observations (c,d) 
during late southern spring in MY 36. The spectral profile along the black line 
in ais shown in e and reveals a marked increase in reflectance up to 20% in the 
BLU filter when frost is present. Errors are from the uncertainty in the absolute 
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calibration of the instrument and are about -3% (ref. 43). The illumination 
direction is indicated by the arrows in the bottom right corner of each image. 
North is up in all panels. The CaSSIS image IDs are shown in order (a-d): 
MY35_008465_192 0_NPB, MY36_020297 350.3 NPB, MY36_020366_190_1 NPB 
and MY36_020478 190 3 NPB. Credit: a, ESA/TGO/CaSSIS under a Creative 
Commons license CC-BY-SA 3.0 IGO. 


pole-facing slopes” *’. However, the presence of frost at the tropics 
(-0° N latitude) was not expected because of higher average surface 
temperatures and lower humidity”. Some studies predicted that 
on most of Mars’s surface, small amounts of H,O frost can condense 
nightly if radiative cooling is strong enough**”’”. For example, extremely 
small amounts of water frost have been observed to condense near the 
equator on the high thermal emissivity calibration targets of NASA’s 
Opportunity rover”. 

Most of the Martian atmosphere is composed of CO, gas, and 
therefore CO, frost can also form if surface temperatures are low 
enough”. On the basis of nightly surface temperatures and thermal 
modelling, it was shown that in the equatorial regions CO, frost may 
condense diurnally”. Predictions’? indicated that putative CO, 
frost deposits may persist for only a few minutes after sunrise (~6:00 
Local Solar Time (LST)) before sublimating back into the atmosphere. 
Follow-up global surveys, utilizing early-morning colour observa- 
tions from the Thermal Emission Imaging System (THEMIS*) were 
conducted to search for these frosts in the equatorial regions, but no 
evidence of morning CO, frost was identified™*. 

Observations by the Colour and Stereo Surface Imaging System 
(CaSSIS*) on board the European Space Agency’s (ESA's) Trace Gas 
Orbiter (TGO) provide strong evidence for morning frost deposition 
on the equatorial Tharsis volcanoes. We present here these obser- 
vations coupled with supporting evidence from other instruments 
and modelling. 


Observations of frost 

Early-morning images (LST = 7:11; latitude = 18.5° N, longitude = 
-133.5° E; spatial resolution = 4.5 m pixel’) of Olympus Mons caldera 
acquired by CaSSIS (at UTC 2022 November 25) in the late northern win- 
ter (solar longitude (Ls) -345°) on Mars year (MY) 36 first revealed bluish 
deposits (at -500 nm) on sections of the caldera floor and rim (Fig. 1). 
The CaSSIS observation suggests a spatial correlation between the blu- 
ish deposits and topography (Fig. 1d). The deposits are concentrated 
on the caldera floor but are absent on well-illuminated warm slopes 
and farther northon the volcano flank. The finding was confirmed five 
days later with a High Resolution Stereo Camera (HRSC)* observation 


acquired on 2022 November 30 (LST = 7:20; latitude = 18.2° N, longi- 
tude = -133.2° E; spatial resolution = 800 m pixel”), which revealed 
that the diffuse bluish ‘halo’ deposit was ubiquitous on the entire cal- 
dera floor and rim (Fig. 1b,c). The halo is absent on the volcano flanks 
and is concentrated only at the mountain summit. During the CaSSIS 
detection, the Nadir and Occultation for Mars Discovery (NOMAD”) 
spectrometer was operating and acquired a ride-along observation 
(instantaneous field of view = 17.5 km x 0.5 km). The nadir spectral data 
acquired inthe NOMAD limb nadir and solar occultation (LNO) channel 
revealed that the deposit is frost (Fig. 1e) as indicated by the elevated ice 
index values (more than 30 confidence; Methods and Supplementary 
Figs. 1and 2). 

Repeat imaging by HRSC shows that the frost deposits on top 
of Olympus Mons (Fig. 1b) appear only in the early Martian morning 
(LST = ~7:00-7:30; latitude = 18.2° N, longitude = —133.2° E) and are 
spatially correlated with a geological bright halo unit (Extended Data 
Fig. 1a,d). This unit may be dust that is relatively brighter than the sur- 
rounding material due to different grain size or texture”. This bright 
halo unit is also observed in Context Camera” images (Extended Data 
Fig. 2a). Materials consisting of smaller particles may exhibit different 
thermophysical properties such as lower thermal conductivity*° and 
high thermal emissivity”. Surfaces with such properties cool down 
moreat night and warm up moreslowly inthe morning, further enhanc- 
ing the likelihood and duration of frost formation. This latter point is 
illustrated by CaSSIS observations of frost on dust deposits that have 
not been removed by winds (Extended Data Fig. 2b,c). As shown by CaS- 
SIS, frost may also condense leeward of small craters where air-fall dust 
can accumulate and is perhaps less compact (Extended Data Fig. 2d). 
Porous and less-compact materials provide more nucleation sites for 
frost formation”. Outside of the bright halo, frost is found near the 
northern rim of Olympus Mons, but its emplacement is more localized 
(Extended Data Fig. 2e-i). In conclusion, the observed frost patterns on 
Olympus Mons, particularly in areas with geologically distinct bright 
dust deposits, underscore the importance of thermophysical proper- 
ties such as low thermal conductivity and high thermal emissivity, as 
well as surface texture, in governing the formation, distribution and 
persistence of frost on Mars. 
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Fig. 3 | Seasonal and diurnal trends of frost on Olympus and Arsia Montes. 
a-d, Rose diagrams showing the seasonal (a,c) and diurnal (b,d) frost detections 
by CaSSIS over Olympus Mons (a,b) and Arsia Mons (c,d). The width of each bin 
reflects the number of CaSSIS observations. Frost is detected around northern 


b H Frost 
E Nofrost 


12h 


6h 


spring equinox (Ls = -0°) on Olympus Mons and around southern winter solstice 
(Ls = ~90°) on Arsia Mons. Frost is detected only in the early-morning hours 
(-7:00-8:00 LST). The negative detections in the early morning bins correspond 
to observations that were acquired in warmer seasons. 


Within the CaSSIS database, 13 instances of frost have been found 
(Extended Data Fig. 3). These include detections not only on the larg- 
est Tharsis volcanoes of Olympus, Ascraeus and Arsia Montes but 
also on the smaller-sized Ceraunius Tholus shield volcano (Extended 
Data Fig. 4). In one case, the frost deposits on Arsia Mons (Fig. 2a) are 
observed in the early Martian morning and during southern winter 
solstice (LST = ~8:00, Ls = ~90°, latitude = -8.7° N, longitude = -121.1° E). 
The frost line dividing warm and shadowed slopes is, however, not 
observed in the repeat CaSSIS observations of this location, which were 
acquired during late southern spring (Fig. 2b—d). Photometric analysis 
shows that frost is associated with an increase in ratioed reflectance of 
up to 20% at wavelengths (Fig. 2e; CaSSIS blue (BLU) filter bandwidth 
is 390-570 nm (ref. 43)). The fact that frosty surfaces are sometimes 
brighter only at blue wavelengths, implying a lower spectral slope, 
can also be observed in a linearly stretched CaSSIS BLU filter image 
(Extended Data Fig. 5 and Methods) and average spectra from a k-means 
clustering analysis**** applied on a topographically corrected and 
photometrically normalized CaSSIS cube (Supplementary Fig. 3). The 
photometric and clustering analyses suggest that the frost deposits 
are probably very thin. 

CaSSIS observations of Olympus and Arsia Montes indicate 
diurnal and possibly seasonal trends in frost deposition (Fig. 3). The 
four detections in Olympus Mons (Fig. 3a,b) are clustered around 
the early-morning hours (LST = ~7:00-7:30) and northern spring 


equinox (Ls = ~-320-40°). Similarly, the four detections in Arsia Mons 
(Fig. 3c—d) fall within a slightly wider time range (LST = ~ 7:00-8:30) but 
around the southern winter solstice (Ls = ~45-145°). The early-morning 
non-detections in Arsia Mons fall within the southern summer period, 
which suggests seasonality (Extended Data Fig. 6b), but the lack of 
early-morning observations in the northern summer precludes us 
from making the same conclusion for the detections in Olympus Mons 
(Extended Data Fig. 6a). We removed observations at extremely high 
solar incidence angles (>85°) because of low image signal-to-noise ratio 
(SNR), and therefore there is an observational bias towards LSTs at about 
6:00 (Methods). Collectively, CaSSIS observations suggest that the 
frost cycle over Martian volcanoes is ephemeral and exhibits variability 
on multiple timescales. It appears to be influenced by diurnal patterns, 
probably reflecting daily temperature fluctuations. In addition, there 
is a probable control by the Martian seasons, indicating a longer-term 
variation in the frost cycle. On the basis of the CaSSIS observations, 
while there are indications of diurnal and seasonal influences on frost 
deposition on Martian volcanoes, these observations alone cannot 
definitively determine the composition of the frost. Therefore, we use 
simulations of surface temperatures as a proxy for frost composition. 


Surface temperatures indicate water frost 
At the time of CaSSIS frost detections in Olympus Mons (Fig. 1) and 
Arsia Mons (Fig. 2), the surface temperatures calculated by the general 
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Fig. 4 | Modelled surface temperatures at Olympus Mons and Arsia Mons. 

a,b, Annual surface temperatures at four different local mean solar times (LMST). 
c,d, Diurnal surface temperatures at Olympus Mons (Ls = 350°) (c) and 

Arsia Mons (Ls = 90°) (d) as predicted by the GCM. Inc,d, blue and red horizontal 
dashed lines depict CO,” and H,O frost point, respectively. On both volcano 
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calderas at the time of CaSSIS image acquisition, CO, frost point is not reached. 
This indicates favourable conditions for H,O ice. The simulations were 


conducted at geographical coordinates 18.75° N, -133.75° E for Olympus Mons 
and -8.75° N, -121.25 °E for Arsia Mons. 


circulation model (GCM*) via the Mars Weather Research and Forecast- 
ing (WRF”) model are inconsistent with CO, frost. The stability of CO, 
frost at higher altitudes necessitates exceptionally low temperatures, 
specifically below 140 K, to maintain its solid state*°. The predicted 
surface temperatures (at -150 km model resolution) are ~150 K and 
185 K at -7:00 LST in Olympus Mons and at ~8:00 LST in Arsia Mons, 
respectively (Fig. 4). In addition, advanced high-resolution mesoscale 
modelling, with a model resolution of 5.47 km, reveals a substantial 
temperature difference between the surface temperature and the local 
CO, frost point at the locations of CaSSIS observations, with a differ- 
ence of approximately 10 K at Olympus Mons (Fig. 5d) and over 55 K at 
Arsia Mons (Extended Data Fig. 7d). In fact, the surface temperatures 
predicted at each CaSSIS frost location (Table 1) consistently exceed 
the CO, frost point, corresponding to the mean surface temperature 
~162 K (excluding C3 and C6). The stratification of water vapour in 
Mars’s atmosphere, especially near the surface, is not well understood**, 
making the determination of the H,O frost point challenging due to 
its considerable variability; however, it is generally accepted that this 
point occurs at around 180 K (ref. 14). Since the predicted surface tem- 
peratures at the time of CaSSIS, HRSC and NOMAD observations are 
too warm, this suggests that CO, frost is unlikely, hence providing sup- 
port for the presence of water frost. At these seasons (Ls = 346.7° for 
Olympus Mons and Ls = 93.8° for Arsia Mons), CO, frost was also not 
observed by the Thermal Emission Imaging System” or by the Emirates 
Mars InfraRed Spectrometer”. Interestingly, the GCM also predicts that 
some CO, frost may be present at Ls = ~O-150° and at around sunrise 
(5:00-6:00 LsT) in Arsia Mons (Fig. 4b). This result is consistent with 
previous studies indicating CO, frost formation from minutes to tens 


of minutes after sunrise in the equatorial regions”. However, such 
potential CO, frost deposits would sublime very quickly and would 
be difficult to detect by cameras and spectrometers due to low SNR™. 
In addition, we investigated the possible role of CO, frost in regolith 
gardening and slope streak formation on Mars*°**“’. We found no slope 
streaks on the calderas of the largest Tharsis volcanoes or any obvious 
differences in talus boulder shapes and sizes (Methods and Extended 
Data Fig. 8). These results suggest that the diurnal CO, or H,O frost 
cycle plays a minor (if any) role in landscape evolution at these sites. 


Microclimate and water ice amount 

Our high-resolution mesoscale simulations reveal the distinct microcli- 
matic conditions induced by the topography of the Tharsis volcanoes, 
as shown in Fig. 5 and Extended Data Fig. 7. Specifically, within the 
calderas of Olympus Mons and Arsia Mons, we observe a substantial 
reduction in surface atmospheric pressure and near-surface horizontal 
wind speeds compared with the surrounding areas. For example, within 
the caldera of Olympus Mons (Fig. 5b), the atmospheric pressure is 
estimated at only 110 Pa, compared with 160 Pa at the mountain’s base. 
Similarly, in the area of Arsia Mons (Extended Data Fig. 7b), the pres- 
sure is about 100 Pa, notably lower than the over 200 Pa found in the 
adjacent plains. Moreover, the near-surface horizontal wind speeds 
within Olympus Mons (Fig. 5c) are estimated at less than10 ms“, in stark 
contrast to the approximately 30 ms! observed along the volcano’s 
flanks. In the case of Arsia Mons (Extended Data Fig. 7c), the wind speeds 
arebelow5 ms ‘within the caldera, compared with roughly 20 m s*on 
the flanks, highlighting the profound impact of volcanic topography 
on localized weather patterns. 
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Fig. 5 | Microclimatic conditions simulated over Olympus Mons. cause noticeable variations in local pressure, wind velocities and temperature 
a-d, Utilizing MarsWRF high-resolution mesoscale modelling (at the time of gradients. The black outline across all panels highlights the boundary of the 
CaSSIS observation Fig. 1d), this figure presents the influence of Olympus Mons’s Olympus Mons caldera, while the black dashed rectangle marks the area 
topography on its local climate, as shown by elevation gradients (a), surface observed by CaSSIS, as referenced in Fig. 1d. The CO, frost point in the area of 
atmospheric pressure (b), near-surface horizontal wind patterns (c) and the CaSSIS observation is exceeded by about 10 K. By contrast, the CO, frost point in 
deviation in temperature between the Martian surface and the local CO; frost Arsia Mons is exceeded by around 60 K (Extended Data Fig. 7d). 
point (d). The topography of the Olympus Mons caldera is demonstrated to 
Table 1| CaSSIS frost detection times and predicted surface temperatures by the GCM 
ID Locality Ls (°) LsT (h) Estimated surface Estimated CO, frost Estimated surface H,O ice stability? 
temperature (K) temperature (K) pressure (Pa) 
Cl Olympus Mons 326.98 7:35 160.83 136.83 114.03 Yes 
C2 Arsia Mons 93.94 8:04 186.83 136.81 113.60 Yes 
C3 Arsia Mons 148.07 7:25 216.70 136.77 112.83 a 
C4 Olympus Mons 338.96 TAU 156.77 136.69 111.34 Yes 
C5 Olympus Mons 35.24 6:58 164.34 136.58 109.35 Yes 
C6 Arsia Mons 43.39 8:35 215.77 137.90 136.63 z 
C7 Ascraeus Mons 35.20 7:06 162.75 138.41 148.78 Yes 
c8 Olympus Mons 344.08 7:12 153.04 136.67 110.98 Yes 
c9 Ceraunius Tholus 16.82 7:05 No data No data No data No data 
C10 Ascraeus Mons 42.80 7:20 166.16 138.31 14619 Yes 
ci Arsia Mons 44.19 7:08 155.19 138.68 155.59 Yes 
C12 Olympus Mons 18.80 6:10 147.07 136.61 109.79 Yes 
C13 Olympus Mons 42.84 7:13 170.53 136.55 108.69 Yes 


*Frost is found only on shadowed slopes. 
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Furthermore, our GCM simulations suggest that the thickness of 
water frost deposits is on the order of 1 um (Methods). However, this 
estimate carries considerable uncertainty due to the unknown quanti- 
ties of water-vapour-column abundances. To refine this estimate, we 
reference radiative transfer calculations”, which suggest a minimum 
thickness of 100 um, while laboratory experiments” imply a thickness 
of about 10 um (Methods). By adopting the median thickness of 10 pm 
for the water frost, and considering that the frost deposits are confined 
tothe calderas of Olympus, Arsia, Ascraeus Montes and Ceraunius Tho- 
lus, we estimate that there is a transfer of approximately 1.5 x 10° kg of 
water ice between the surface and the atmosphere (Methods). 


Possible sources of water vapour 

The seasonal trends as shown by the set (n = 13) of CaSSIS observations 
suggests an atmospheric phenomenon driven by water transport due 
to large-scale seasonal changes, such as sublimation of the seasonal 
ice cap in the opposite hemisphere and transportation of humid air 
into the volcano calderas by upslope winds. Seasonal processes have 
been observed at a wide range of Martian latitudes” and may also 
apply to the Tharsis region. For example, the activity of the Aphelion 
Cloud Belt peaks at Ls ~40°-140° (ref. 54), and in general little cloud 
activity is observed at Ls -245°-320°. Similarly, afternoon orographic 
clouds have been detected by Mars Reconnaissance Orbiter’s Mars 
Color Imager” over the Tharsis volcanoes”. The seasonal observation 
of water-vapour enrichment over Tharsis” shows increased abundances 
around northern spring equinox (Ls 0°), consistent with the CaSSIS 
detections of frost close to this season in Olympus Mons. Therefore, we 
hypothesize that this water-vapour enrichment” may be the source of 
the frost deposits detected in our study. The transport of water vapour 
from high latitudes to the Tharsis highlands could be facilitated by 
large-scale atmospheric eddies”. This process could be further aug- 
mented by strong upslope winds, driven by a combination of thermal 
effects and mountain gravity waves”, facilitating the movement of 
moisture over the volcano calderas. The local topography-induced 
circulation” and microclimatic conditions within the caldera (shown 
in Fig. 5 and Extended Data Fig. 7) may create favourable conditions 
for water frost condensation during the cold Martian nights. Within 
these calderas, ~150,000 tons of water ice is exchanged daily between 
the regolith and the atmosphere during the cold Martian seasons. 
Although this amount is relatively a small fraction of the seasonal inven- 
tory of water vapour in the Martian atmosphere (-10” kg) (ref. 14), 
it is important in the context of localized Martian environmental pro- 
cesses. Understanding these micro-environments is crucial for acom- 
prehensive understanding of Mars’s hydrological cycle. 

Itis conceivable that dormant volcanoes can emit CO,, water vapour 
and minor amounts of SO, (ref. 58) via diffuse outgassing from the rego- 
lith?’ If the observed water frost deposits are of volcanic origin, their 
distribution may constrain models for present-day outgassing fromthe 
interior. However, on Mars, SO, has not been detected” and no thermal 
hotspots have been found”. A volcanic source for the condensate can- 
not completely be ruled out, but further tests for trace species (CO,, 
H-S and SO,) would be useful to explore the likelihood of this potential 
mechanism. Consequently, we conclude here that the newly detected 
frosts on Tharsis volcano calderas are probably of atmospheric origin. 
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Methods 

CaSSIS frost observations 

We surveyed ~4,200 CaSSIS images (acquired up to 2022 February 05) 
with illumination geometries of 50-90° incidence within dusty, low 
thermal inertia (<100 TIU) regions (60° N—30° S). Only images that 
include the latest CaSSIS radiometric and absolute calibration were 
used in this study**>*. 

The images used in this study consisted of early (6:00-9:00 LST) 
and late (15:00-18:00 LST) times. Analysis and comparison in these two 
local time regimes may help the distinction between early-morning 
and late-afternoon phenomena. During the survey, it was noticed 
that most CaSSIS images acquired at extremely high solar incidence 
angles of 85-90° contain colour and calibration artefacts due to the 
decrease in SNR and/or an increase in aerosol contribution from the 
atmosphere’. Consequently, the images with colour artefacts were 
labelled as ambiguous and were not used for further analysis. 

Frost detections relied on the use of CaSSIS NPB (near infrared 
(NIR) = 940, panchromatic (PAN) = 670, blue (BLU) = 497 nm) and 
synthetic RGB (red-green-blue; PAN and BLU only) products. These 
filter configurations allow a convenient separation between frosty 
and frost-free terrains. In CaSSIS colour products, frosty areas appear 
bluish, and/or whitish, and sometimes are bright only in the BLU filter 
(relative to frost-free areas; also see Supplementary Figs. 9 and 10). In 
support, we observe bluish frost deposits in HRSC colour images shown 
in Fig. 1b and Extended Data Fig. 1 (composites of blue (440 nm), green 
(530 nm) and red (750 nm) channels). 

As shown by previous studies” deposits are usually correlated 
with topography (prefer poleward-sloping terrains). Therefore, if 
both conditions were met (colour and topographic correlation), it was 
considered a strong indication of surface frost. As a final procedure, 
each of these candidate detections was then analysed using a spectral 
profile tool in the Environment for Visualizing Images software. This 
procedure extracts the pixel irradiance over flux (//F) values between 
two manually selected points crossing the potentially frosted regionin 
each filter. The profiles were then normalized by a mean //F of a nearby 
frost-free, relatively flat region of interest (ROI), a well-established 
method to cancel out some of the atmospheric and topographic 
effects*”°°”°. If the frost deposits were brighter in the BLU filter than 
the surrounding frost-free terrains by at least 3% (within CaSSIS abso- 
lute uncertainty*’), then such images were flagged as potential frost 
detections. This survey yielded many frosty sites (not shown here) at 
latitudes -40° N and ~30° S. However, because these latitude bands are 
dominated by known seasonal frost deposits”’”*” and we do not have 
a robust method to distinguish between seasonal and diurnal frost, 
we further narrowed our filtering criteria. The final frost detections 
analysed here were restricted to equatorial -20° N to ~10° S latitudes 
(outside of the seasonal mid-latitude regions). In this work, only equa- 
torial sites that included visible evidence of frost are considered. 

The spectral profile shown in Fig. 2e was computed by dividing 
each pixel along the profile by an average pixel value extracted from an 
ROlin Extended Data Fig. 5d. The ROI (>100 pixels in size) was selected 
ona frost-free and relatively flat terrain as suggested by the low slope 
values in the CaSSIS digital elevation model of this site. CaSSIS digital 
elevation models were produced by a pipeline developed at the Astro- 
nomical Observatory of Padova, National Institute for Astrophysics ™”. 


NOMAD-LNO spectral processing 

The NOMAD instrument is a suite of three high-resolution spectrom- 
eters also on board TGO, offering nadir infrared observations through 
its LNO channel”. This channel covers the 2.2-3.8 um spectral range 
where several spectral features of ice are distributed over different 
wavelengths. Nevertheless, the NOMAD-LNO spectrometer has the 
particularity of not observing the entire spectral range at once. The data 
are acquired through small spectral windows, representing specific 
diffraction orders of the diffraction grating. Each LNO observation can 


select amaximum number of 6 diffraction orders every 15 seconds to 
ensure the best possible SNR”. The LNO footprint (instantaneous 
field of view) is 17.5 km x 0.5 km (ref. 75), which provides enough spa- 
tial scale to resolve the caldera of Olympus Mons. In this work, we use 
spectrally and radiometrically calibrated LNO data converted into a 
reflectance factor. The 2.7 um ice band is the strongest inthe LNO spec- 
tral range, resulting from both CO, and H,O ice absorption. Although 
the use of this band is not suitable for quantifying the amount of ice 
(easily saturated), it is effective for detecting homogeneous deposits 
(both CO, and H,O ice), as demonstrated with the ice index value”. 
This spectral parameter uses two diffraction orders. It is based on the 
combination of high reflectivity at continuum wavelengths with a more 
pronounced absorption in the 2.7 um band. Initially defined as the 
spectral ratio between the reflectance factors of order 190 (continuum 
part, 2.32-2.34 um) and 169 (short wavelengths shoulder of the 2.7 um 
band, 2.61-2.63 um) (ref. 77), we adjust the ice index by considering 
the available orders of the joint CaSSIS-NOMAD observations, that 
is, orders 190 and 168 (2.64-2.65 um). 

In nadir mode, the variability in the reflectance factors is caused 
mainly by the surface albedo variations resulting from the different 
absorption of the Martian surface mineralogy’ ®. To remove spatial 
albedo variations over the explored Martian surface, we normalize the 
LNO reflectance factors to the Martian albedo. The adjusted ice index 
(II) can thus be defined as: 


l= R233 i R263 
OMEGAz2;2 | OMEGA 62 


where R,is the LNO reflectance factor value averaged around the central 
wavelength iof the LNO spectrum, fitted by a third-degree polynomial 
to mitigate the spectral oscillations resulting from the instrumen- 
tal characteristics of the LNO channel, which become significant on 
the edges of each order (Supplementary Figs. 1 and 2). OMEGA, is the 
OMEGA albedo map” based on reflectance spectra in the near infrared 
as NOMAD-LNO R;. Two OMEGA albedo maps are used in this work: one 
defined at 2.32 um for order 190 and the other defined at 2.62 um for 
order 168. Studies have shown that this spectral parameter identifies 
spatially extensive and abundant ice deposits when the index values 
are three sigma higher than their average value over ice-free mid- 
latitude terrain’. 


Mars GCM modelling 

We perform the Martian GCM simulation for the entire MY 36 using the 
MarsWRF model, whichis the Mars adaptation of the general-purpose 
planetary atmosphere model, planetWRF”. Here the GCM set-up is 
based ona previous study*ć examining the Martian planetary boundary 
and dust-turbulence interaction over a decade, from MY 24 to MY 34, 
which hosted three global dust storms. The reference model set-up*° 
was validated against NASA’s Mars Climate Sounder (MCS) observa- 
tions on board the Mars Reconnaissance Orbiter, radio occultation 
observations from ESA’s Mars Express orbiter, as well as the in situ 
observations from NASA’s Mars Science Laboratory Curiosity rover. 
This model set-up consists of a semi-interactive two-moment dust 
transport model” within the MarsWRF framework, in a way that the 
dust is lifted, mixed by model winds and sedimented, as guided by 
observed maps of column-integrated dust optical thickness**’. Via 
this method, model processes govern the vertical dust distribution 
and related dust radiative heating, yet the horizontal dust distribution 
is guided to match the orbiter observations. In this model, the hori- 
zontal dust distribution is constrained to follow observations. In this 
model, the two-stream correlated k-distribution scheme is used for the 
short-wave and long-wave radiative transfer**. We use a Mars-specific 
boundary-layer turbulence parameterization scheme, which allows 
us to obtain the surface-atmosphere exchange coefficients® Surface 
properties of the MarsWRF model, such as the topography, albedo, 
emissivity and thermal inertia, are acquired from the datasets of the 
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Mars Orbiter Laser Altimeter’? and Thermal Emission Spectrometer 
(TES”*) observations, where the details are presented in another study”. 
Here we increased the horizontal model grid spacing ofthe GCM from 
5° x 5° to 2.5° x 2.5° (ref. 46), enabling better spatial coverage to provide 
more realistic boundary and initial conditions to our mesoscale simula- 
tions. We used 52 vertical sigma layers extending up to the model top 
of 100 km. The predicted surface temperatures are shown in Table 1. 

Our modelling methodology is based on a previous study by 
MarsWRF*°°°, Mesoscale simulations for Fig. 5 and Extended Data 
Fig. 7 were forced with initial and boundary conditions acquired by 
GCM simulations corresponding to the same seasonal conditions of 
CaSSIS observations shown in Figs. 1 and 2. The plots we present in 
terms of winds, pressure and temperature correspond to the local 
hours of observations. We nested three mesoscale domains in our GCM 
domain (see Supplementary Fig. 12 for details). Mesoscale domains use 
prescribed boundary conditions, derived either from GCM predictions 
(as in the case of d2) or from another mesoscale domain (d3 and d4). 
The GCM grid has a horizontal resolution of approximately 150 km. 
We progressively increased the horizontal resolution with a factor of 
three for our nested mesoscale domains. Our innermost domain, d4, 
has a horizontal resolution of 5.47 km. To assess the accuracy of our 
mesoscale predictions, we compared MarsWRF surface temperature 
predictions with the surface temperature observations by MCS and TES 
available for Olympus Mons and Arsia Mons regions at around 3:00 LST 
(Supplementary Fig. 13). We considered a sufficient Ls range of MCS 
and TES observations (Ls 310-360 for Olympus Mons and Ls 75-100 
for Arsia Mons) to provide a sufficient set of observations to acquire a 
temperature map to be compared with MarsWRF simulations. These 
observations range from 1:00 LST to 4:00 LST, and MarsWRF estima- 
tions at the corresponding local times are compared for validation. 
The modelled surface temperatures for Olympus Mons caldera are 
within 10 K of the observations and within a few degrees Kelvin for 
Arsia Mons. It isimportant to note that these predictions carry uncer- 
tainties, particularly in regions with complex topography such as the 
Tharsis volcanoes. 


Surface frost thickness and mass estimations 
The MarsWRF GCM incorporates the phase transition and transport 
mechanisms of water vapour and ice, facilitating a parameterization 
of the Martian hydrological cycle that aligns with the methodologies 
outlined by previous studies*”. This parameterization enables the 
model to approximate the surface frost layer thickness to about 1 um at 
the locations in our study. However, it is important to acknowledge the 
inherent uncertainties associated with such estimations, particularly 
due to the limitations of physical parameterizations within Martian 
atmospheric models. These uncertainties are most pronounced inthe 
prediction of atmospheric variables in regions lacking empirical obser- 
vational data, suchas the deposition rates of atmospheric volatiles. 

In a recent experimental investigation, one study” systemati- 
cally evaluated the interaction between water frost deposition and 
the optical properties of a Martian soil simulant, specifically Mars 
Global Simulant (MGS-1°°). The experimental design involved the 
controlled deposition of water frost on the surface of the simulant, 
followed by precise measurements of both the spectral reflectance 
and the thickness of the frost layer. The findings indicate that a frost 
layer thickness ranging from approximately 10 to 20 um is required 
to significantly attenuate the characteristic red slope of the spectral 
reflectance, aligning with the observed morning frost brightening in 
the blue wavelengths by approximately 10-20% as detected by the 
CaSSIS instrument. Furthermore, the study demonstrates that a rela- 
tively thin frost layer of about 100 pm is sufficient to flatten the visible 
spectrum, effectively neutralizing the spectral features. 

Radiative transfer models” can provide an additional constraint 
on the frost thickness estimation via the minimum optical depth (T) 
necessary for frost visibility at CaSSIS visible and LNO near-infrared 


wavelengths. For example, witha tof 107, weanticipatea minimalimpact 
on albedo, less than 0.1 at CaSSIS visible wavelengths and negligible at 
LNO near-infrared wavelengths, given the single-scattering co-albedo 
is around 10~ for visible light and less than 10“ at 2.6 um. However, LNO 
observations indicate a discernible albedo reduction at near-infrared 
wavelengths, suggesting a higher optical depth than 10”. This implies 
that the frost’s grain radius and/or thickness must exceed 5 pmand1 um, 
respectively. Ifthe grain radius is about 1 um, then the frost layer’s thick- 
ness could be significantly greater, approximately 100 um. 

To conduct a preliminary quantification of the frost mass, 
we assumed a uniform frost layer thickness across all identified 
frost-covered regions, as observed by CaSSIS. The geographical extent 
of the frost coverage was approximated to the combined surface areas 
of the calderas of Martian volcanoes such as Arsia Mons, Olympus 
Mons, Ascraeus Mons and Ceraunius Tholus. By integrating the uni- 
form frost thickness with the delineated area and adopting the density 
value for pure ice, we derived an initial estimate of the total frost mass. 
This approach provides a rudimentary yet insightful approximation 
of the frost mass, acknowledging the broad-scale estimative nature 
of this calculation. 


Boulder size measurements 

To investigate a potential effect of the diurnal frost cycle on the overall 
geomorphology and landscape evolution, we studied the shape of 
mass-wasted boulders across six sites of interest. Here we compare the 
sizes of boulders on volcanoes with frost as determined by CaSSIS (two 
sites in Olympus Mons and one in Arsia Mons) and on volcanoes where 
frost has not been detected (Tharsis Mons, Jovis Tholus and Ulysses 
Tholus). Because frost accumulates preferentially on poleward-facing 
slopes on Mars”, here we focused only onnorth-facing and south-facing 
slopes. This might reveal whether there are considerable differences 
in boulder sizes due to frost weathering”. 

We used eight map-projected High Resolution Imaging Science 
Experiment (HiRISE)*° images in Geographic Information System 
(QGIS) to determine the three principal dimensions of each identified 
boulder. The first dimension is defined as the longest distance between 
two points on the boulder as visible from orbit. Similarly, the second 
dimension is defined as the diameter of the boulder orthogonal tothe 
first dimension. Last, the third dimension is defined as the height of the 
boulder as estimated using shadow length and solar incidence angle. 
Intotal, we identified and measured 63 boulders across the six sites. All 
derived measurements were plotted on ternary diagrams” using the 
Tri-Plot software”. These diagrams relate the three principal dimen- 
sions of each boulder, visualizing its overall shape as well as similarities 
and differences within and across the studied sites. 


Data availability 

CaSSIS data can be found on the University of Bern repository 
(https://observations.cassis.unibe.ch/) and the ESA’s Planetary Science 
Archive (https://archives.esac.esa.int/psa). NOMAD-LNO observations 
are also found on the ESA’s Planetary Science Archive. 


Code availability 
The PlanetWRF model for Martian GCM and mesoscale simulations is 
accessible by request at https://planetwrf.com/. 
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Frost halo 


5.7) 

Extended Data Fig. 1| Diurnal variations of frost halo on Olympus Mons. deposit is likely composed of fine-grained dust with low thermal conductivity, 
HRSC images of Olympus Mons (lat =18.2°N, lon = -133.2°E) acquired at different which facilitates frost formation. North is up in all panels. HRSC image 

local times in MY 36. (a, b) Late morning images showing no evidence of frost on IDs: hn705_0000 (a), hn772_0000 (b), hn889_0000 (c), hn948_0000 (d). 

the bright halo deposit surrounding the volcano caldera. (c, d) Early morning Credit: ESA/DLR/FU Berlin. 

images revealing the presence of frost on the bright halo deposit. The bright halo 
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Extended Data Fig. 2 | Irregular frost distribution on the outskirts of 
Olympus Mons caldera. (a) Bright halo visible in the CTX global mosaic 

(Dickson et al., 2018). The bright halo deposit is also visible in HRSC non- 
detections in Extended Data Fig. 1. (b) Dark windstreak and triangular bright 
unremoved dust deposits seen in CTX. (c) CaSSIS morning observation of 
froststreaks (lat/lon =18.28°N, -134.24°E) that correlate with bright dust deposits 
seen in CTX. (d, e) Frost deposits leeward of small craters (lat/lon: 18.27°N, 
-134.24°E and 18.83°N, -133.65°E respectively). (f, g) Froststreaks that are parallel 
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Frost trap 


Froststreak 


to local winds (lat/lon: 18.90°N, -133.65°E and 18.9°N -133.74°E respectively). 
(h, i) Small frost deposits on the rim of a collapse pit (lat/lon: 18.95°N, -133.69°E) 
and on the levee of a lava channel (lat/lon: 19.00°N, -133.71°E). North is up in 

all panels. CaSSIS IDs: MY37_023825_162_0 (LST =7:13AM; Ls = 42.84; c,d) and 
MY36_015229 _160_0 (LST = 6:57AM; Ls = 35.24 ; e-i). Credit: a,b, NASA/JPL/ 
MSSS/The Murray Lab; c-l, ESA/TGO/CaSSIS under a Creative Commons license 
CC-BY-SA 3.0 IGO. 
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Extended Data Fig. 3 | CaSSIS frost detections in the Tharsis volcanic region. Tholus volcano. Frost has not been observed yet on Pavonis Mons and other 
Tharsis volcanoes. The basemap is the color hillshade MOLA data at 64 pixels per 


Frost was detected only on and around the calderas of the three largest volcanoes 
such as Olympus, Arsia and Ascraeus Montes, but also on the smaller Ceraunius degree resolution. Credit: NASA/JPL/GSFC. 
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Extended Data Fig. 4 | Frost on the caldera floor of Ceraunius Tholus volcano. different local time featuring no frost. Both CaSSIS images in (b) and 
(a) Wide angle view of Ceraunius Tholus (lat = 24.0°N, lon = —97.1°E) with CaSSIS (d) are acquired at similar incidence (and phase) angles, which suggests 
early morning observation overlain on the CTX mosaic. (b) Zoomed in view of (a). that photometric effects are not the cause of surface blueing. North is up 


White rectangle marks the close up in (c). (c) Ubiquitous frost coverage on the in all panels. CaSSIS image IDs in order: MY37_023134_024_3_NPB (a-c) and 
caldera floor and the apparent absence of frost on the caldera rim. MY36_022599_024_0_NPB (d).a,b,d, NASA/JPL/MSSS/The Murray Lab; 
(d) CaSSIS color NPB image of the Ceraunius Tholus caldera acquired ata a, b, c, d, ESA/TGO/CaSSIS under a Creative Commons license CC-BY-SA 3.0 IGO. 
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Spectral profile 
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eee! DEM 
fr = i sca 
Extended Data Fig. 5 | Greyscale CaSSIS filter images of a small crater (Methods) overlaid. The frost deposits are brighter and visible in the CaSSIS BLU 
presented in Fig. 2. (a-c) Three individual filters: NIR (860-1100 nm), PAN filter. North is up in all panels. CaSSIS ID: MY35_008465_192_0 (a-c) and DEM ID: 
(550-800 nm) and BLU (390-570 nm). (d) CaSSIS DEM of the same scene with the CAS-DTM-MY36_020366_190_1-OPD-03-01 (d). Credit: ESA/TGO/CaSSIS under a 
location of the spectral profile and the frost-free ROI used in the reflectance ratio Creative Commons license CC-BY-SA 3.0 IGO. 
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Extended Data Fig. 6 | Seasonal and local time coverage of CaSSIS 

early morning observations of Olympus Mons and Arsia Mons. CaSSIS 
image coverage over Olympus Mons (a) and Arsia Mons (b). Based on these 
observations frost is not detected during late morning hours in Olympus Mons 


and around southern summer solstice (Ls ~270°) in Arsia Mons. The shaded grey 
region shows the CaSSIS observational bias (from 90 - 85° solar incidence) due to 
the low signal-to-noise ratio. Most observations were discarded in this region due 
to spectral ambiguity. The black lines mark the local sunrise time. 
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Extended Data Fig. 7 | Microclimatic conditions simulated over Arsia Mons. and the local CO, frost point. The complex caldera topography of Arsia Mons is 
This figure illustrates the impact of Arsia Mons’ topography on localized shown to substantially influence local pressure distributions, wind velocities, and 
atmospheric conditions (at 8AM, Ls= 90°, Fig. 2a), depicted through (a)elevation thermal gradients. Notably, surface temperatures at sites identified by CaSSIS for 
gradients, (b) surface atmospheric pressure, (c) near-surface horizontal wind frost presence exceed the CO, frost point by approximately 60 K, suggesting the 
patterns, and (d) the temperature differential between the Martian surface predominance of H,O ice in these frost deposits. 
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Extended Data Fig. 8 | Triangular (ternary) diagrams of boulder shape across 
six sites of interest. Boulder shape and size analysis for Olympus Mons 

(a,b), Arsia Mons (c), Tharsis Tholus (d), Jovis Tholus (e), and Ulysses Tholus (f). 
The slope aspect of each site is indicated on the respective panel. The Olympus 
Mons (a,b) and Arsia Mons (c) sites were found to feature a distinct early 
morning frost signature (positives, blue), the other three sites not (negatives, red) 


(d-f) - yet there is no obvious difference in boulder shape across those sites 


(the colored polygons underline the distribution of points). C = compact, 
P=platy, B = bladed, E = elongated, V = very; a = longest boulder dimension, 
b =intermediate boulder dimension, c = smallest boulder dimension. 
HiRISE image IDs (a-f): ESP_014275_1990_RED, ESP_043272_1980_RED, 
ESP_047439_1990_RED, PSP_009884_1980_RED, ESP_057843 1715 RED, 
ESP_012612 1940_RED, ESP_033711_1985 RED, and ESP_045619 1835 RED. 
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